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ABSTRACT: There are many impedance-based methods to an- 
alyze the oscillation of grid-connected converter system (GCS) 
via dq frame, sequence frame and polar frame, which can be 
transformed into each other and are equivalent in judging the 
stability of the system. However, the discrepancy and adaptation 
of these methods has not been discussed theoretically. Therefore, 
the impedance models of GCS are established in above three 
kinds of coordinates respectively, then the mathematical and 
physical characteristics of them are analyzed. Secondly, the de- 
coupling of impedance models in different frames are analyzed, 
and the differences of the impedance criteria in stability analysis 
and mechanism interpretation are given. Furthermore, the sub- 
synchronous oscillation caused by the phase-locked loop consid- 
ering converter four quadrant operation scenario are illustrated, 
and the existence of open-loop transfer function right half plane 
poles and its sensitivity to uncertain parameters are analyzed 
from the nominal performance and robustness of the impedance 
criteria. On this basis, the adaptation of several impedance crite- 
ria of GCS are analyzed. Theoretical analysis and simulations 
confirm the validity of the proposed conclusion. 


KEY WORDS: impedance-based stability criterion; renewable 
energy; power factor; robustness; nominal performance; 


摘要 : 基于 dq 坐标 、 序 坐标 和 极 坐 标的 阻抗 法 是 研究 变 流 
器 振荡 问题 的 常用 方法 , 不 同 阻抗 法 可 相互 转化 并 在 判断 系 
统 稳定 性 方面 等 价 , 然而 它们 的 差异 性 和 适应 性 尚 待 深入 探 
讨 。 为 此 , 首先 在 上 述 三 种 坐标 系 分 别 建 立 并 网 系统 阻抗 模 
型 ， 并 对 比分 析 其 数学 和 物理 特征 ; 其 次 ， 分 析 了 不 同 坐 标 
系 下 阻抗 模型 的 可 解 耦 性 ,给 出 了 所 导出 的 阻抗 判 据 在 主导 
环节 分 析 和 机 理解 释 方面 的 具体 区 别 。 进 一 步 ， 以 锁 相 环 高 
参与 度 的 次 / 超 同 步 振荡 问题 为 例 ， 针 对 变 流 器 四 象限 运行 
时 的 不 同 场景 , 分 别 从 阻抗 判 据 的 标 称 性 和 重 棒 性 角度 分 析 
了 开 环 传递 函数 右 半 平 面 极 点 的 存在 性 和 对 不 确定 性 参数 
的 灵敏 性 , 并 在 此 基础 上 分 析 了 几 种 现 有 阻抗 判 据 的 适用 场 
景 和 条 件 。 最 后 ， 利 用 仿真 验证 所 提 结 论 的 有 效 性 。 
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0 引言 


目前 大 部 分 的 电力 电子 并 网 装备 通过 锁 相 环 
(phase locked loop，PLD) 与 电网 实现 同步 ， 但 随 着 
新 能 源 装 备 的 大 规模 渗透 ， 交 流 电 网 逐渐 变 弱 ， 新 
能 源 并 网 系统 的 小 扰动 稳定 性 问题 日 益 突出 ， 且 引 
起 的 次 / 超 同步 振荡 逐渐 成 为 制约 新 能 源 发 展 的 关 
键 因素 之 一 中。2009 年 美国 德州 某 风 场 发 生 的 次 
同步 谐振 造成 多 台风 机 援 棒 保 护 电 路 损坏 外、2011 
年 我 国 新 疆 哈 密 发 生 风 电机 群 参与 的 严重 次 同步 
振荡 事故 造成 特 高 压 直 流 功 率 又 降 的 严重 后 果 W]， 
事后 的 研究 表明 ，PLL 是 造成 问题 的 主要 原因 。 

由 于 新 能 源 设备 模型 阶 数 高 、 且 设备 控制 结构 
参数 往往 属于 厂家 机 密 ， 详 细 建 模 难 以 实现 ， 导 致 
基于 时 域 的 状态 空间 法 在 分 析 电 压 源 型 变 流 器 

(voltage-sourced converter, VSC) 并 网 系统 的 小 扰 
动 稳定 问题 时 存在 局 限 。 目 前 工程 上 往往 采用 基于 
频 域 的 阻抗 法 ， 通 过 VSC 和 电网 的 端口 外 特性 进 
行 稳定 性 分 析 ， 其 阻抗 模型 不 仅 易 于 测量 ， 且 阻抗 
建 模 可 将 系统 划分 为 网 络 侧 和 设备 人 出， 适合 解决 
VSC 与 弱电 网 的 耦合 问题 ， 故 受到 了 广泛 关注 区 9]。 

阻抗 建 模 依 赖 于 选取 的 坐标 系 ， 常 见 的 有 基于 

同步 旋转 坐标 系 的 dq 坐标 阻抗 模型 4020 和 极 坐 标 
阻抗 模型 0422， 以 及 静止 坐标 系 下 的 序 阻 抗 模 型 0 
1 。 不 同 坐 标 系 下 的 阻抗 矩阵 均 含 有 非 对 角 耦 合 项 ， 
为 了 得 到 可 以 简单 应 用 的 阻抗 判 据 ， 目 前 主要 的 方 
法 有 两 种 ， 一 种 是 基于 小 增益 稳定 准则 导出 多 输入 
多 输出 系统 的 阻抗 判 据 ， 如 基于 禁止 域 或 范 数 的 判 
据 ， 但 这 类 方法 存在 一 定 的 保守 性 05]。 另 一 种 是 通 
过 数学 变换 或 物理 等 值 构造 单 输 入 单 输出 (single- 
input single-output, SISO) 模型 , 进而 导出 阻抗 判 据 。 
例如 ， 当 VSC 和 负载 均 运 行 在 单位 功率 因数 工 况 ， 
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VSC 侧 的 dq 阻抗 矩阵 具有 主 对 角 占 优 的 性 质 ， 可 
以 名 略 非 对 角 耦 合 项 进行 稳定 性 分 析 09。 

然而 ， 一 些 新 的 并 网 标准 对 VSC 提出 了 无 功 
需求 07, 因此 需 考虑 VSC 功率 多 象限 输出 运行 (后 
文 简称 多 象限 运行 ) 的 情况 19。 此 时 忽略 看 合 项 会 
导致 较 大 误差 ,需要 采用 较 复杂 的 广义 奈奈 斯 特 判 
据 (generalized Nyquist criterion, GNC) 进行 分 析 
0920。 为 了 解决 阻抗 矩阵 耦合 导致 的 稳定 分 析 复杂 
或 存在 误差 的 问题 , 文献 [21] 考 虑 dg 阻抗 模型 的 非 
ARAL, HES TAR q 轴 电 流 的 稳定 判 据 来 分 
析 系 统 的 稳定 性 。 文献 [22] 给 出 了 dq 阻抗 模型 和 序 
阻抗 模型 的 转化 关系 ， 并 考虑 正 负 序 分 量 间 耦合 的 
影响 导出 了 修正 的 序 阻抗 判 据 。 文 献 [23] 提 出 了 将 
dq 阻抗 模型 转化 到 of ER SISO 阻抗 模型 的 方法 ， 
实现 了 阻抗 模型 在 dq 域 与 of 域 形 式 上 的 统一 。 文 
献 [24] 考 虑 功率 因数 变化 ， 导 出 了 基于 极 坐标 系 的 
广义 阻抗 判 据 , 将 VSC 并 网 系统 解 耘 为 SISO 系统 
分 析 。 文 献 [25] 探 讨 了 考虑 外 环 动态 后 基于 广义 阻 
抗 的 系统 稳定 性 分 析 方 法 。 

上 述 研究 均 可 以 将 复杂 系统 的 振荡 问题 简化 
A SISO 系统 的 问题 ， 为 了 厘清 这 些 方法 的 区 别 ， 
己 有 研究 表明 不 同 稳定 判 据 的 适用 范围 和 场景 存 
在 差异 。 针 对 广义 阻抗 判 据 和 序 阻抗 判 据 , 文献 [26] 
从 复 电 路 角度 阐明 了 它们 的 导出 机 制 和 振荡 物理 
解释 ， 证 明了 两 种 判 据 表征 的 稳定 裕 度 在 鲁 棒 性 方 
面 存 在 差异 ， 文 献 [27] 提 出 了 稳定 判 据 的 适用 性 评 
估 方 法 , 阐述 了 稳定 判 据 的 误 用 可 能 会 导致 其 在 所 
分 析 频 段 内 的 稳定 裕 度 失效 和 重 棒 性 不 足 的 问 
题 。 然 而 , 上 述 研究 均 局 限 在 额定 工 况 ， 当 考虑 VSC 
多 象限 运行 时 ， 阻 抗 矩 阵 耦 合 加 剧 将 导致 阻抗 法 的 
稳定 分 析 更 为 复杂 ， 目 前 尚未 有 体系 化 的 量化 评估 
制 分 析 它 们 的 适用 
为 此 ， 本 文 将 文献 [27] 中 的 稳定 判 据 适用 性 评 
估 方 法 拓展 到 VSC 多 象限 运行 的 场景 ， 探 讨 dg AA 
标 、 正 负 序 坐标 和 极 坐 标 阻抗 判 据 的 特点 和 适用 
PE. 首先, 基于 三 种 坐标 系 分 别 构建 VSC 并 网 系统 
阻抗 模型 并 导出 相应 阻抗 判 据 。 在 此 基础 上 ， 通 过 
分 析 不 同 阻抗 判 据 的 耦合 特性 说 明了 它们 在 主导 
环节 分 析 和 机 理解 释 方面 的 具体 区 别 。 其 次 ， 以 锁 
相 环 主导 的 VSC 振荡 问题 为 例 ， 分 别 从 开 环 传递 
函数 的 标 称 性 和 和 鲁 棒 性 角度 对 比 了 不 同 阻抗 判 据 
的 适用 性 。 最 后 ,在 MATLAB /Simulink 中 搭建 VSC 
并 网 系统 模型 ， 通 过 仿真 验证 了 上 述 结论 的 有 效 
性 。 研 究 表明 ， 针 对 VSC 多 象限 运行 场景 下 PLL 
主导 的 小 扰动 稳定 问题 ， 基 于 广义 阻抗 判 据 进行 稳 
定 分 析 和 控制 设计 适用 性 更 好 。 


ie 


1 变 流 器 多 象限 运行 
用 性 问题 
1.1 ” 变 流 器 并 网 系统 阻抗 模型 的 构建 

考虑 图 1 所 示 由 VSC 和 交流 电网 组 成 的 并 网 
系统 ， 其 中 VSC 采用 基于 PLL 的 双环 矢量 控制 策 
略 ， 考 虑 PLL， 外 环 ， 电 流 内 环 和 电压 前 馈 控制 环 
动态 。 图 中 Ew Vabe 和 Wine BY ATCA A HAP. JF 
网 公共 连接 点 和 VSC 输出 端口 三 相 电 压 ，Iabe 和 
Ieabe 分 别 为 流 过 滤波 电感 和 线路 电感 的 三 相 电 流 。 


景 下 阻抗 判 据 的 适 


Ler.Cr 和 工分 别 为 滤波 电感 ,滤波 电容 和 线路 电感 ， 
下 标 “4 ”、“g”* 表 示 电 气量 在 dq 旋转 坐标 系 下 的 d 


轴 和 4 轴 分 量 。 如 不 特别 说 明 ， 文 中 的 dq 坐标 为 
全 局 同步 坐标 ， 并 假设 VSC 和 电网 的 导 纳 和 矩 阵 都 
不 存在 不 稳定 极点 。 此 外 ， 为 简化 问题 ， 本 文 分 析 
过 程 中 首先 考虑 功率 外 环 控制 , VSC 外 环 采 用 不 同 
有 
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1 VSC 控制 框图 
Fig.1 Diagram of the grid-connected converter 


VSC 并 网 系统 的 dq 坐标 导 纳 模型 及 对 应 的 特 
征 方程 为 (具体 变量 推导 详 见 附录 A) 


=| Y, H A - 3 Kaloo AV, (1) 
à 40 do +Y Lao AV, 
Co E 


式 中 
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_ Gy (s) 
1+G(s)G (s) © 1+G, 
Y = -Gru (s)G (s) , r-l; H 

2 1+G(s)G (s)V, Q s 
Yaq we 和 Yaq net 分 别 为 VSC 和 电网 的 dq FN FE BE, 
考虑 功率 外 环 时 认为 GAp= GRp=GPp。 
VSC 并 网 系统 的 序 坐 标 导 纳 模型 以 及 对 应 的 
特征 方程 可 以 表示 为 


Yy 


Yr kr 

det [Yno Yn J= de a a (6) 
式 中 : Yp=[(stiwo)Le]'+Ci(stjwo), YN= [(s-jwo)Ls]'+ 
Ci(s-jwo), Yen 为 序 坐标 下 的 导 纳 矩阵 。 
极 坐 标 系 下 的 导 纳 模型 , 可 以 反映 电压 幅 值 / 相 
角 和 电流 幅 值 / 相 角 间 的 相互 作用 ， 且 和 dq 坐标 导 

纳 模型 间 可 以 相互 转化 I， 满足 如 下 关系 
Ya vse = PY,, ve (7) 
Yo p = PY (8) 


pmj 


Ss 


Rp P= Ce ARE E 


-sino coso 
压 的 相 角 ，YG 为 极 坐标 下 的 导 纳 矩阵 。 
VSC 并 网 系统 的 极 坐标 导 纳 模型 以 及 对 应 的 
特征 方程 可 以 表示 为 


Al Yo 0 AV 
=| P+ (9) 
IAO, 0 -¥,+%, |J VAO 
和 


Settee a 
TAO, L, EAG, | |VA@, 


F (11) 
Yi Y, 
式 中 : 91 和 Ov 为 极 坐标 下 电流 和 电压 矢量 的 角度 。 
1.2” 变 流 器 并 网 系统 阻抗 模型 的 结构 特征 

当 VSC 在 多 象限 运行 时 ， 极 坐标 系 下 VSC 并 
网 系统 的 特征 方程 为 


Yy 0 l -1 
det! I, +Y,,P+P| —F!+C,F |$=0 (12) 
027.47, L 


det cao + YG net | = s 


结合 Yvi Yv2 和 Yv; RR E HAG oO) 2 所 
AN» E PLL 参与 度 高 的 中 低频 段 ，Yvs 和 Yv; 的 幅 
值 远大 于 Yvyi。 此 时 Yo vse 为 主 对 角 占 优 矩 阵 , 说 明 
极 坐标 系 下 设备 侧 导 纳 模型 在 数学 上 可 近似 解 耦 ， 
物理 特点 为 设备 的 电气 量 幅 值 和 相 角 相互 解 耦 。 
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Fig.2 Bode diagram of the Yvi, Yv2 and Yv3 
此 外 ， 将 (12) 中 设备 侧 不 含 PLL 动态 的 传递 函 
数 系 卫 折 算 到 网 络 侧 ， 则 设备 侧 动态 将 严格 解 耦 ， 


而 网 络 侧 矩阵 转化 为 如 下 形式 


= 
OO 
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Yam =P[ Ket LE +c | (13) 
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式 (13) 满 足 和 矩阵 主 对 角 元 素 相 同 ， 非 对 角 元 素 
相反 的 结构 ， 因 此 极 坐标 系 下 VSC 并 网 系统 的 特 
征 方程 可 以 表示 为 ， 


(| 
det + =0 (14) 
0 Y| |b a 


式 (14) 中 设备 侧 为 对 角 阵 , 网 络 侧 为 对 称 阵 。 这 
种 特殊 结构 有 利于 解决 传统 阻抗 分 析 方 法 在 中 低 
频段 由 于 阻抗 矩阵 的 非 对 角 耦 合 项 导致 的 分 析 复 
杂 或 存在 较 大 误差 的 问题 ， 可 以 简化 稳定 分 析 [133。 

对 于 dq 坐标 阻抗 模型 , 结合 (1)-(3) 可 知 当 VSC 
满足 = Att, 基于 dq 坐标 和 极 坐 标 系 的 特征 方 
程 等 价 。 然 而 ,考虑 VSC 多 象限 运行 的 场景 (140#0)， 
VSC 和 电网 的 阻抗 矩阵 均 存 在 非 对 角 耦 合 项 , 其 设 
备 侧 和 网 络 侧 导 纳 和 矩阵 均 含 有 非 对 角 元 耦合 项 。 

而 对 于 式 (4)-(5) 所 示 的 序 阻抗 模型 ， 虽 然 网 络 
侧 正 序 和 负 序 阻抗 间 可 以 实现 解 耦 , 但 PLL 的 存在 
导致 设备 侧 正 负 序 阻抗 在 中 低频 段 仍 存在 较 大 的 
耦合 项 ， 如 果 忽 略 耦 合 项 可 能 会 导致 稳定 分 析 结 果 
错误 L493。 此 外 ， 当 考虑 振荡 频率 较 高 的 稳定 问题 
如 线路 LC 谐振 29， 由 图 2 可 知 wi 的 幅 值 远大 于 
Yva 和 Yv3, 此 时 序 阻抗 模型 设备 侧 和 网 络 侧 导 纳 算 
阵 可 近似 对 角 解 耦 ， 而 dq 坐标 和 极 坐 标 阻抗 模型 
中 设备 侧 Yui 的 影响 不 可 和 忽略。 


1.3 ”多 种 阻抗 模型 的 适用 性 问题 


本 文 探讨 基于 dq 坐标 、 序 坐标 和 极 坐标 阻抗 
判 据 在 分 析 VSC 接 入 弱 网 后 次 / 超 同步 振荡 问题 时 


4 中 国 电 


的 适用 性 。 三 种 阻抗 模型 可 以 相互 转化 且 对 应 的 特 
征 值 相 同 (证明 详 见 式 (29))， 所 以 在 稳定 性 判别 结 

果 上 一 致 。 然 而， 结合 1.2 节 的 分 析 ， 考 虑 VSC 多 
象限 运行 的 场景 ， 三 种 坐标 系 下 阻抗 模型 结构 存在 
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A yA 
el 卷 
a 


点 a 的 周 数 。 
2) dq 坐标 阻抗 判 据 
文献 [21] 通 过 选择 并 网 电流 的 9 轴 分 量 得 到 等 
效 SISO AZ, 并 基于 VSC 的 9 轴 阻 抗 和 等 效 9 $ 


差异 ， 物 理 上 反映 了 三 种 模型 的 输入 输出 电气 量 耦 
合 特性 不 同 ， 这 导致 基于 模型 导出 的 阻抗 判 据 在 稳 
定 机 理 分 析 和 裕 度 量化 方面 存在 区 别 〈 详 见 第 四 章 
算 例 ), 但 目前 对 此 尚未 有 合理 的 评估 方法 , 导致 在 
有 具体 场景 下 难以 合理 选择 相 适 应 的 阻抗 分 析 方 法 。 
由 此 引出 如 下 问题 : 为 了 对 PLL 高 参与 度 次 / 
超 同 步 振荡 问题 进行 有 效 分 析 和 控制 ， 阻 抗 判 据 不 
仅 应 该 反映 系统 是 否 稳定 ， 也 需要 准确 量化 系统 的 
稳定 裕 度 ， 揭 示 控 制 动 态 对 稳定 性 的 影响 规律 。 考 
E VSC 多 象限 运行 场景 ， 基 于 dq 坐标 、 序 坐标 和 
标的 阻抗 判 据 在 分 析 PLL 与 电网 的 交互 作用 、 
量化 系统 稳定 裕 度 方面 的 适用 性 应 该 如 何 评估 ? 
接 下 来 分 别 基于 三 种 阻抗 模型 导出 阻抗 判 据 ， 
从 频 域 控制 角度 探讨 它们 的 特点 和 适用 性 。 


2 阻抗 判 据 的 特征 分 析 


2.1 典型 阻抗 判 据 的 回顾 

VSC 多 象限 运行 场景 下 ， 基 于 1.2 节 的 分 析 可 
知 dq 阻抗 窍 阵 d, q 轴 回 路 相互 耦合 且 不 可 忽略 ， 
目前 常见 的 分 析 方 法 有 两 种 : 1) 对 系统 的 回 比 矩阵 
进行 相似 对 角 变 换 , 导出 基于 特征 轨迹 分 解 的 GNC 
判 据 中。2) 保留 d 轴 或 g 轴 回 路 ， 通 过 和 矩阵 的 舒 
尔 补 变换 将 其 余 物 理 量 折算 消去 来 构建 等 效 的 
SISO 模型 并 导出 阻抗 判 据 R11。 

序 坐 标 系 下 可 以 通过 修正 序 阻抗 模型 分 析 系 
统 的 稳定 性 23， 极 坐标 系 下 VSC 多 象限 运行 时 阻 
抗 矩 阵 满足 式 (14) 所 示 特 殊 结构 ， 可 通过 广义 阻抗 
的 等 价 变换 降 阶 特征 方程 导出 广义 阻抗 判 据 031。 

1) GNC 判 据 

dg 坐标 系 下 基于 GNC 判 据 分 析 ， 首 先 需 要 建 
WAZA IB] EKER 


> 


u 


LaS Vid Ve (15) 
通过 相似 对 角 变换 得 到 Lag 的 特征 函数 为 
A(s)=HL,,J (16) 


式 中 : A(s)=diag[Ai(s), A2(s)], H 和 J 分 别 为 La 的 
左右 特征 向 量 和 矩阵 。 

GNC 判 据 : 若 Lag 存 在 p 个 右 半 平 面 (right half plane, 
RHP) 极 点 ， 当 其 特征 轨迹 A1(s), 42(s) 在 复 平 面 逆 时 
针 包 围 (-1y0) 的 总 周 数 为 p 时 闭环 系统 稳定 ， 记 作 


Veneta, =p (17) 
式 中 ,enc[Z,a] 代 表 闭 合 曲线 二 逆 时 针 包围 复 平 面 一 


m 


电网 阻抗 的 传递 函数 分 析 系 统 稳定 性 。 
当 保 留 系统 的 q 轴 变 量 , 式 (3) 可 以 转化 为 如 下 
特征 方程 
det[ Yi ]det| 到 eae (yey) "Ero (18) 

4 Yi, 18) 
det| x (人 ve =0 (19) 


即 原 系统 的 稳定 问题 转化 为 消去 d 轴 的 输入 / 
输出 变量 得 到 等 效 9 轴 阻 抗 模 型 的 稳定 性 问题 ， 进 
而 得 到 如 下 特征 方程 


=I 
Yer + Yq Ya" ) Yt =0 (20) 
SF 


Xam 


AP, YaN Ye TIA q HH VSC 和 网 络 侧 的 导 纳 ， 
且 满 足 ， Ye+yY 。 
当 满 足 det[Yqqm]z0 时 ， 得 到 基于 等 
抗 的 稳定 判 据 Cda 坐标 阻抗 判 据 ) 
lat, =0 (21) 
式 中 : L5 Ya q [Yagm]! ,代表 qË VSC 阻抗 和 等 效 q 
轴 电 网 阻抗 的 阻抗 比 。 
dq 坐标 阻抗 判 据 : AL, FE p 个 RHP 极点 ， 则 当 
其 频率 特性 曲线 在 复 平 面 逆 时 针 包 围 (-1, j0) 的 周 数 
为 p 时 闭环 系统 稳定 ， 记 作 
enc| L,,-1|=p (22) 
3) 序 阻抗 判 据 
基于 正 负 序 坐 标 系 的 序 阻抗 模型 ， 可 以 利用 修 
正 序 阻抗 判 据 分 析 系 统 的 稳定 性 [3, 选择 保留 正 序 
、 JE AVP, PFE A 内 折算 消去 ， 可 以 构 
等 效 的 SISO 模型 并 导出 修正 序 阻抗 判 据 ( 负 序 
后 文 简称 为 序 阻抗 判 据 : 
(er. 20 (23) 
式 中 , Le=YpvscZPnet, YPnet=Ype?”,Zpnet=|Ypnet|!, Ypvse=Yvi 
eF? to Yvat {[loYvot(YvitYuyje] + Lo ¥v3/2-loYv2y!t 1, 
Ly=YyvscZNnet,YNnet=Yne’”,Znnet=| YNnet] !, Ynvsc=Yvie® + 
ToYv2t{ [loYv2t+(Yv1t+YP )e7?] li+[JoYv3/2-1o Yv2] yte 
序 阻抗 判 据 : 若 Lb 存在 p 个 RHP 极点 ,LN 存在 m 
个 RHP 极点 , 则 当 Lb 和 LN 均 满足 奈 硅 斯 特 稳定 判 
据 ， 则 并 网 系统 稳定 。 记 作 
enc|L,,—-1]=p 
< [5,-1]=m A 
4) 广义 阻抗 判 据 
基于 极 坐标 系 的 广义 阻抗 模型 ， 可 以 利用 广义 


效 q AARE 


阻抗 判 据 分 析 稳 定性 [49。 联 立 式 (9)-(11)， 令 
[AE,EA0E]'=[0,0]'〈 认 为 无 穷 大 电网 稳定 )， 并 将 与 
功率 因数 相关 的 传递 函数 矩阵 YviP 折算 到 网 络 侧 ， 
得 到 如 下 特征 方程 


Y, 0 AV 
D YP +Y, asa |=) (25 
| 0 | 0 -Y,, ex. Pl vi + G_net J eA ( ) 


选择 保留 式 (25) 中 的 相 角 变量 Abov, 将 幅 值 变 
te AV 折算 消去 ， 可 以 构建 等 效 的 SISO 模型 并 导 
出 广义 阻抗 判 据 


1+L, =0 (26) 
sh, Lo=YoavscZanet, Yovsc= lo(Yv3-Yv2)/2, Zonet= 
[Yonet]!, Yoner= -oYvo/2 + { [loYv2 + (Yvı+Yp)e” J! 

+ [loYvot(Yvit Yo J”! y! 
广义 阻抗 判 据 : 若 Lo i CARA ae Fe, JS 
网 系统 稳定 。 由 于 Lo 不 存在 RHP 极点 ， 当 Lo 的 频 
率 特性 曲线 不 包围 (-1,0) 时 , 系统 小 扰动 稳定 。 记 作 
enc[L,,—1]=0 (27) 


2.2 ”阻抗 判 据 的 耦合 特性 分 析 

1) GNC 判 据 

由 式 (28) 可 以 看 出 基于 特征 轨迹 的 GNC 判 据 
是 将 原始 的 2 维系 统 解 丰 为 两 个 虚拟 单 变量 系统 。 
如 图 3 Ca) 所 示 。 原 系统 的 动态 响应 是 两 个 虚拟 单 
变量 系统 的 组 合 , 其 加 权 系数 由 Laa 左右 特征 向 量 
的 乘积 决定 ， 有 具体 为 


A 0 
1+4 
0 Ay 
1+4 
由 于 不 同 坐 标 系 可 通过 线性 变换 相互 转化 ， 因 

此 基于 不 同 阻抗 模型 求解 的 特征 值 相同 。 以 dq AK 
标 和 极 坐 标 阻抗 模型 为 例 给 出 证 明 如 下 : 

0= det| Ly, ~ 2al | 

=det| P'L,P-P'(A,,1)P | 

= det| P~ |det| Le — Al |det[P] (29) 

= det| Le - 2,1] =0 

=> Ao = Aaa 
式 中 :46 和 444 分 别 为 极 坐标 和 dq AA tr BESTE REKI 
特 征 值 , Lo=YG vsc/YG net 为 极 坐标 阻抗 回 比 矩 阵 o 
因此 ， 同 一 个 系统 在 不 同 的 坐标 系 下 建 模 ， 再 
基于 特征 轨迹 进行 稳定 分 析 , 结论 完全 一 致 。 然 而 ， 
车 于 特征 轨迹 的 稳定 分 析 方 法 难以 得 到 直观 的 特 
征 函数 显示 表达 式 ， 系 统 参数 和 稳定 性 的 关系 较 难 
刻画 ， 不 利于 指导 参数 设计 05。 同 时 ， 基 于 特征 轨 
迹 的 控制 设计 需要 考虑 每 文 虚拟 单 变量 系统 的 动 
态 响应 ， 在 工程 上 应 用 较为 复杂 R231。 


L,,[I+L.,] =H" J” (28) 
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2) dq 坐标 阻抗 判 据 

在 单位 功率 因数 工 况 下 基于 dq 坐标 和 极 坐 标 
系 的 阻抗 模型 相同 ， 此 时 dq 坐标 阻抗 判 据 和 广义 
阻抗 判 据 等 价 。 然 而 考虑 VSC 多 象限 运行 场景 
Gaoxt0)，d 坐标 阻抗 判 据 所 分 析 的 开 环 传递 函数 
中 YY 和 Yam AE PLL 的 动态 , 为 便于 说 明 这 种 
分 布 对 稳定 分 析 的 影响 ， 忽 略 外 环 和 前 馈 的 低 通 滤 
BATA, dq 坐标 阻抗 判 据 所 分 析 的 等 效 SISO 
系统 如 图 3 (b) 所 示 ， 发 现 PLL 的 闭环 传递 函数 
GPLL(s) 与 电网 动态 及 设备 运行 的 功率 因数 相互 看 
合 ， 导 致 对 开 环 传递 函数 进行 稳定 分 析 时 较 复 杂 。 
3) 序 阻 抗 判 据 
基于 序 阻抗 的 SISO 系统 如 图 3 (c) 所 示 ， 序 阻 
抗 模型 将 VSC 的 负 序 回路 动态 通过 舒 尔 补 变换 折 
算 到 正 序 回路 , 可 以 发 现在 PLL 参与 度 较 高 的 中 低 
频段 ，PLL 的 闭环 动态 仍然 难以 与 电网 动态 及 设备 
运行 的 功率 因数 相互 解 耦 分 析 。 此 外 ， 基 于 序 阻 抗 
判 据 的 分 析 方 法 需要 考虑 w 从 -oo 到 +oo 变 化 时 的 频 
率 特性 P9， 从 而 增加 了 分 析 的 工作 量 。 

4) 广义 阻抗 判 据 
基于 广义 阻抗 的 SISO 系统 如 图 3 (AMR, 
义 阻 抗 判 据 不 含 功率 因数 引起 的 阻抗 耦合 项 ， 系 统 
的 小 扰动 稳定 问题 可 表示 为 VSC 和 网 络 相 角 的 交 
互 作用 。 进 一 步 结 合式 (26) 可 知 ， 广 义 阻抗 判 据 所 
分 析 的 开 环 传递 函数 中 与 VSC 运行 点 相关 的 电气 
量 均 可 以 显 式 表达 ,PLL 的 全 部 动态 只 出 现在 与 功 
率 因数 无 关 的 Yovsc F, M Zoner 包含 了 系统 其 余 
环节 的 动态 。 因 此 ， 基 于 广义 阻抗 判 据 可 以 通过 
Yovsc 和 Zoner 的 伯 德 (Bode) 图 来 直观 的 分 析 不 同 运 
人 


从 4 
9 1 维 


y 1+ Gru (5)Gi(s) wo Llo sing | sf 1 | 
| sL, Gp, (8)G,(s)J, cos o 
\ J te 


(d) 广义 阻抗 判 据 
3 等 效 SISO 系统 的 数学 模型 
Fig.3 Diagram of the equivalent SISO system 
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3 ”阻抗 判 据 的 适用 性 分 析 


3.1 阻抗 判 据 的 标 称 性 分 析 

当 阻 抗 判 据 所 分 析 的 开 环 传递 函数 工 存 在 RHP 
极点 时 ， 在 稳定 性 分 析 和 镇 定 系统 时 会 受到 诸多 限 
制 ， 可 能 难以 保证 系统 的 标 称 性 要 求 P9]。 
对 广义 阻抗 所 分 析 的 开 环 传递 函数 Lo, 因 VSC 
自身 是 稳定 的 , 则 Lo YovscZaner 不 存在 RHP 极点 
等 价 为 Zoner MIFE RHP 极点 。 

记 ZaNer=NMG)/D(G)， 原 问题 转化 为 Zoner 的 分 
D(s) 没 有 RHP 零点 。 考 虑 到 在 次 / 超 同 步 频段 内 
环 与 前 馈 滤 波 的 传递 函数 满足 Gve(s)~0, Gi(s)~ 
1251, t D(s) 可 以 简化 为 


D(s)=2+4Y,[e (s+ jm)+e™” (s- jm)|L, (30) 
将 Ws 具体 表达 式 代 入 式 (30) 可 以 得 到 
D(s)= As? +Bs+C (31) 


TRH, A=(IoL eKpcosg+T pa), B=(1+ JoLeKicosp t+mK?), 
C=mKi, m=1- woL losing, Kp 和 KK 为 外 环 控制 参数 。 
由 静态 稳定 极限 可 知 jwoLglosing| 小 于 1, 并 且 在 


% dg 阻抗 La x 序 阻抗 le % 广义 阻抗 Le 
s,=-7.792+ j93.96 


I,g=-0.90 10=0.44 


工程 学 R 


A Se. 
第 AE 
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典型 参数 下 IsKpcoso| <|Teo| 恒 成 立 ， 因 此 式 (31) 
多 项 式 系数 A,B,C 总 为 正 ， 根 据 赫 尔 维 获 判 据 B9， 
式 (30) 不 含 RHP 零点 。 因 此 ,广义 阻抗 判 据 对 应 的 
开 环 传递 函数 Lo 不 存在 RHP 极点 。 

然而 ， 在 PLL 参与 度 较 高 的 次 / 超 同步 振荡 频 
Be, 如 果 考 虑 VSC 多 象限 运行 , 较 强 的 频率 耦合 作 
用 可 能 会 导致 dg 坐标 和 序 坐 标 系 下 构建 的 开 环 传 
BERA Ly 和 Le 存在 RHP 极点 , 不 利于 系统 的 稳定 
分 析 。 通 过 如 下 案例 说 明 此 问题 ， 其 中 fprr=15Hz 
CKPprrp=3 rad/s,Kpiu= 9000 rad/s 习 ， 系 统 其 余 参数 如 
附录 B 案例 1 所 示 。 考 虑 VSC 四 象限 运行 场景 

(VSC AY d 轴 、9 轴 输 出 电流 分 别 为 ro= 士 0.90， 

7ao= 士 0.44， 系 统 的 主导 振荡 模式 为 91), WEN Las 
Lp 和 Lo 的 极点 分 布 和 对 应 的 奈 奎 斯 特 曲 线 如 图 4 
所 示 , 可 以 发 现 基 于 三 种 阻抗 判 据 所 对 应 的 RHP 极 


点 数目 和 频率 特性 曲线 均 可 以 准确 判断 系统 的 稳 
定性 , (EEF L Le E VSC 多 象限 运行 的 场景 下 


可 能 存在 RHP 极点 ,导致 系统 稳定 裕 度 难以 量化 ， 
不 利于 指导 控制 设计 。 


IE lq 


em lp wm LO 


To=0.90 To=0.44 $=-1.659+/j99.17 


75 a sy a xx x 
È ola RHP 极 点 Ë We 
轴 
Bee sax ky 
-150 ® 
10 -5 0 7 i ; 
实 轴 实 轴 
150 
15 Bra x 
虚 0 RHP 极 点 虚 RHP 极 点 
轴 T BT Oe Ness le T 
“|. “| ee xix | x 
ao et I peg S 


实 轴 实 
I, =-0.90 To=-0.44 — s,=-1.067+j85.10 


实 轴 实 轴 
70=0.90 To=-0.44 —s,=4.590+ 89.43 


4 变 流 器 四 象限 运行 下 三 种 判 据 的 开 环 极点 分 析 


Fig.4 RHP pole analysis ofthree impedance criteria considering four quadrant operation of VSC 


3.2 ”阻抗 判 据 的 鲁 棒 性 分 析 

阻抗 判 据 不 仅 可 以 判断 系统 是 否 稳定 ， 也 能 通 
过 量化 系统 稳定 裕 度 来 指导 控制 参数 整定 。 传 递 函 
数 回路 整形 是 一 种 常用 的 控制 设计 方法 ， 通 过 校正 
阻抗 判 据 所 分 析 的 开 环 传递 函数 工 的 频率 特性 达到 
期 望 的 控制 目的 ， 为 了 避免 系统 稳定 裕 度 在 受到 不 
确定 性 摄 动 后 急剧 变化 ,L 应 具备 一 定 的 鲁 棒 性 P29。 

系统 的 稳定 裕 度 可 以 通过 工 的 频率 特性 曲线 与 
(-1,.70) 的 最 短 距 离 来 量化 站， 具体 为 

My =1+L(ja,) 

RP: wi 是 主导 振荡 模式 s1 的 频率 。 


(32) 


受到 扰动 人 后 系统 稳定 裕 度 变化 量 可 表示 为 

AM, =|L(ja,)-L(j(@+A))| (33) 

可 以 利用 在 相同 扰动 下 的 AMy 评估 开 环 传递 

函数 的 鲁 棒 性 ,AMv 更 小 的 传递 函数 对 不 确定 性 摄 

动 的 放大 作用 更 小 ， 其 鲁 棒 性 更 好 。 为 描述 开 环 传 

递 函 数 对 不 确定 参数 的 灵敏 度 ， 引 入 环 路 增益 灵敏 
度 指标 的 模 值 27 

L(j@,) 有 L'(s,) 


lo, (s, ) m Lia, j zs)| L'(s, J 


AP: soto: 表示 系统 的 主导 振荡 模式 ， 上 标 
“ ”代表 求 导 运算 ， 有 具体 参考 文献 [27]. 


(34) 


第 期 KIIR e EREKET: 锁 相 环 主导 下 变 流 器 振荡 问题 的 阻抗 分 析 方 法 : 标 称 性 与 鲁 棒 性 7 


式 (34) 可 以 反映 开 环 传递 函数 工 对 系统 振荡 模 进一步 采用 3.2 节 的 案例 分 析 考 虑 不 同 外 环 控 
式 的 灵敏 度 ， 该 指标 越 小 ， 说 明 L 的 鲁 棒 性 越 好 ， 制 下 阻抗 判 据 的 鲁 棒 性 ， 表 2 是 Lg. Lp 和 Lo 分 别 
即 在 考虑 不 确定 性 摄 动 后 其 稳定 裕 度 的 变化 越 小 。 E VSC 四 象限 运行 Ua=0.80, 40.60) 时 的 
步 ， 利 用 式 (34) 所 提 的 量化 指标 评价 La en a ae VSC 

hy fe be, ike H J 义 阻 抗 回 路 的 灵敏 度 最 小 ,所 表征 的 鲁 棒 性 更 好 。 
Pa asia na elie Racers 表 ?直流 电压 外 环 控制 下 阻抗 判 据 的 鲁 棒 性 分 析 
2 YAW AR SLE BX o FE F> 


7 B Tab. 2 Robustness analysis considering DC voltage control 
象限 输出 运行 时 La Lp 和 Lo 的 |OL(s1)| 指 标 值 , 可 以 PAEAS FFM 
RIL VSC 多 象限 运行 场景 下 广义 阻抗 回路 的 灵敏 Taos oso | oo| owa | oa f om 


a pea | om | 060 | ome | os | 0288 
Bia, MEIEN EPEAT asas] oso | oso | ons | osi | oza 
He 1 新 能 源 四 象限 运行 场景 下 阻抗 判 据 的 鲁 棒 性 分 析 -5.702475.75i | -0.80 | -0.60 | 0.0836 | 0.247 | 0.207 
Tab. 1 Robustness analysis of different stability criteria 
网 和 fr K 
RARR POM [dq E] FP 基于 以 上 分 析 ，VSC 外 环 采用 不 同 有 功 /无 功 


2,660 485.63] 0.90 | -044 | om | 0330 | 0.122 空 制 模式 组 合 , 在 分 析 PLL 高 参与 度 的 次 超 同步 振 
| DSO | ee eee 荡 时 不 同 阻抗 判 据 的 标 称 性 和 和 鲁 棒 性 结论 仍然 适 


-16.26+90.13i | -0.70 | | 0.71 | | 0.0669 | | 0.295 | 0.121 $ 7 
i | -0.60 | -0.80 | 0.0660 | 0.312 0.0889 用 ° 本 文 侧重 于 分 析 由 PLL 与 电网 强 耦 合 导 致 的 次 


xd ste ts SN ` / 超 司 步 振 荡 问 题 ， 限于 篇 幅 原 因 ， ` 同 外 环 控制 模 

3.3 NEINA Eh N l E H! 
不 同 外 环 a I Een AERAN 性 分 析 武 组 合 下 的 具体 分 析 不 再 详细 展开 。 此 外 ， 针 对 外 
研究 PLL 高 参 度 的 振荡 异 式 时 ， 避 保留 号 环 和 锁 相 环 参与 程度 较 高 的 多 控制 环节 交织 失 稳 


PLL 状态 强 相关 的 变量 而 将 其 余 次 要 变量 消去 ”。 问题 ， 不 同 阻抗 判 据 的 标 称 性 和 鲁 棒 性 将 在 未 来 的 
由 于 外 环 动态 对 该 振荡 模式 的 参与 程度 较 小 ， 因 此 工作 中 进一步 研究 。 

当 VSC 外 环 采 用 不 同 有 功 /无 功 控制 模式 组 合 时 ， 弱 网 下 PLL 高 参与 度 的 次 / 超 同 步 振荡 问题 ， 
不 同 阻抗 判 据 的 适用 性 分 析 结 果 与 3.1、3.2 节 保 持 可 以 解释 为 相 角 等 效 阻尼 不 足 导 致 vsC 输 出 电流 


一 致 。 结合 如 下 案例 说 明 该 问题 ， VSC 外 环 采 用 定 相 角 振荡 或 失去 同步 ， 其 主导 输出 变量 为 电流 相 
直流 电压 控制 Kpcpr1.2, Koc™l0), VSC 运行 功率 。 条 pm。 因此， 基于 保留 电流 相 角 动 态 构建 的 广义 


因数 为 0.9 w090, -04 ARERI] 阻抗 可 以 更 好 的 反映 系统 的 同步 失 稳 机 理 ， 在 
节 案 例 一 致 。 系 统 主导 特征 值 为 s 一 3.58+j87.9,， W VSC 多 象限 运行 等 复杂 场景 下 可 以 准确 表征 系 
明 系 统 小 扰动 不 稳定 , La LA LMR eaa O ET AEEA 
率 特性 曲线 如 图 5 所 示 ， 其 中 La LA Laaa a O 

有 2、1 和 0 个 RHP 极点 ， 所 对 应 的 奈 奎 斯 特 曲线 ee ee 


稳定 性 判别 上 等 价 。 当 功率 因数 为 1 时 ， 广 义 阻抗 
和 dq 坐标 阻抗 判 据 的 标 称 性 和 和 鲁 棒 性 等 价 。 当 考 
È VSC 多 象限 运行 场景 ,广义 阻抗 判 据 在 分 析 PLL 
参与 度 高 的 次 / 超 同 步 振荡 问题 时 具有 PLL 动态 和 


分 别 包 围 (-1，j0) 点 0 圈 、1 圈 和 2 Fl, 说 明 三 种 阻 
抗 判 据 的 稳定 性 判断 结果 准确 ， 但 由 于 Lo。 和 Lp FF 
在 RHP 极点 ， 导 致 系统 的 稳定 裕 度 量化 困难 。 


x Lq * PP X Lo — J) m Jp ww [0 
5 


150 系统 其 余 环节 动态 相互 解 得 、 所 分 析 的 开 环 传递 函 
15 数 不 存 在 RHP 极点 和 和 鲁 棒 性 更 好 等 优点 。 然而 , 广 
ki 义 阻 抗 判 据 的 适用 场合 也 具有 局 限 性 。 例 如 ， 该 方 
B 法 在 针对 由 电容 主导 的 高 频 振荡 和 线路 LC 谐振 等 

电气 谐振 问题 时 适用 性 不 足 ， 而 序 阻抗 判 据 在 分 析 
ns 这 类 问题 时 表现 出 更 好 的 适用 性 2s29 。 因 此 明确 广 


义 阻抗 判 据 的 适应 场合 至 关 重 要 ， 运行 
5 直流 电压 外 环 控制 下 阻抗 判 据 的 开 环 极点 分 析 ane 判 据 的 适应 场 VSU FARE 
Fig.5 RHP pole analysis considering DC voltage control Ast PA EARRAK 3 所 示 。 

表 3 变 流 器 多 象限 运行 场景 下 不 同 阻 抗 分 析 方 法 对 比 


Tab.3 Comparison with the different impedance-based method considering multi quadrant operation of VSC 


广义 阻抗 分 析 方 法 序 阻抗 分 析 方 法 dg 阻抗 分 析 方 法 
全 局 坐标 系 对 应 的 极 坐标 正 负 序 坐标 全 局 坐标 系 


GNC 判 据 


广义 阻抗 判 据 序 阻抗 判 据 dq 坐标 阻抗 判 据 


稳定 判 所 enc|L,,-l]}=p 2 
定 enc| L,(s),—-1]=0 ie enc| L (s),-1]=p Danks) =D 
振荡 物理 意义 VSC 和 电网 相 角 阻抗 谐振 | VSC 和 电网 序 电路 谐振 | VSC 和 等 效 q 轴 电 网 阻抗 谐振 / 


8 中 国 电 


在 振荡 频段 可 忽略 
内 环 前 馈 滤 波 环节 动态 


解 耦 条 件 


函数 性 质 高 频段 


TIR 
鱼 棱 性 


4 算 例 分 析 与 验证 
4.1 稳定 性 分 析 与 验证 

由 于 广义 阻抗 判 据 能 直观 描述 VSC 多 象限 运 
行 场景 下 PLL 与 系统 其 余 环 节 的 交互 作用 , 所 以 接 
下 来 据 此 研究 此 场景 下 PLL 动态 与 系统 稳定 性 的 
关系 。PLL 的 闭环 传递 函数 Grrr 为 典型 二 阶 系 统 
35， 首先 讨论 不 同 PLL 的 闭环 带宽 fer 和 阻尼 比 ¢ 
对 VSC 并 网 系统 稳定 性 的 影响 。VSC 运行 功率 因 
BA 0.8 (140=0.80，Joo= 一 0.60),， 系统 参数 如 附录 B 
所 示 。 图 6(a) 给 出 三 种 ferr 下 (保持 0.1) Yovsc 
和 Yoner ÑY Bode 图 ,其 fett=15Hz (KprL1p=18 rad/s, 
Kpr11=8100 rad/s”); forr=20Hz (KpLip=24 rad/s, 
Keri 14400 rad/s’); fert=25Hz (KPLLP=30 rad/s, 
Kpii=22400 rad/s”); 可 以 发 现 PLL 闭环 带宽 主要 影 
响 Yovsc 谐振 尖峰 的 频率 ， 随 着 fori 增加 ，Yeovsc 与 
Yoner 的 交点 频率 增加 ， 且 相位 裕 度 降低 。 图 6(b) 
给 出 了 不 同 5 下 (hz 保持 15Hz) Yovsc 和 Yoner 的 
Bode 图 ， 其 中 C=0.1 (KPLLP=18 rad/s, KpriL1=8100 
rad/s’); ¢=0.15 (KpL1ip=27 rad/s, Kpi11= 8100 rad/s’); 
(=0.2 (Kptip=36 rad/s, KpL11=8100 rad/s’); 可 以 看 出 


of Ni A 


(a) PLL 带宽 对 稳定 性 的 影响 


m= COSP=1 mm COSP=0.8 m= COSP=0.6 


(c) 功率 因数 对 稳定 性 的 影响 


标 称 性 中 频段 不 存在 中 频段 可 能 存在 
环 传递 | RHP 极点 RHP 极点 


鲁 棒 性 较 好 


机 工程 学 


在 振荡 频段 可 忽略 
锁 相 环 和 外 环 动态 
VSC 多 象限 运 
可 能 存在 RHP 极点 
单位 功率 因数 下 


报 


单位 功率 因 玫 
运行 工 况 下 可 


行 场景 下 


频段 鲁 棒 性 较 


好 


解析 复杂 ， 难 以 


观 表现 


PLL 阻尼 比 主 要 影响 Yovsc 谐振 尖峰 的 峰值 ， 相 同 
fo F PLL 的 阻尼 比 越 小 ， 谐 振 尖 峰 越 大 。 


其 次 ,研究 功率 因数 和 


的 影响 。 


Yonet 的 Bode 


电网 强度 对 系统 稳定 性 


图 6(c) 24H 


8 了 不 


同 功 率 因 


数 下 Yovsc 和 


数 为 1 Qal, Iyo=0), 


Yovsc 和 Yoner 的 频率 特性 


曲线 不 相交 , 系统 的 稳定 


mmm, =15Hz =a /11=20Hz em f=25Hz 


裕 度 较 大 。 当 功率 因数 为 0.8 (71ao=0.80, fuo= 一 0.60) 
和 0.6 时 〈7uo=0.60， Iqo=—0.80), Yavsc 和 Yoner 的 
幅 值 曲线 相交 ， 且 相位 差分 别 为 167* 和 178°, RAR 
稳定 裕 度 逐渐 减 小 ,存在 一 定 的 振荡 风险 。 图 6(d) 
给 出 不 同 电 网 强度 下 Yovsc 和 Yoner 的 Bode Kl, “4 
线路 电感 由 0.4pu 增加 至 0.67pu 的 过 程 中 ， 随 着 电 
网 强度 下 降 , 系统 的 稳定 裕 度 减 小 , 失 稳 风险 增加 。 
假设 在 Els 时 VSC 端 电压 受到 扰动 ， 在 上 述 
参数 下 VSC 的 输出 电压 响应 曲线 如 图 7 所 示 ， 系 
统 的 稳定 性 与 频 域 分 析 结 果 相 吻合 。 以 上 分 析 说 明 
不 同 的 参数 变化 可 能 会 导致 弱 网 下 PLL 与 系统 其 
余 环节 交互 作用 加 剧 , 在 处 理 PLL 高 参与 度 的 小 扰 
动 稳定 问题 时 , 利用 广义 阻抗 判 据 能 够 使 PLL 动态 
与 其 余 环节 相互 解 耦 , 便于 直观 分 析 PLL 与 系统 的 
交互 作用 并 揭示 失 稳 机 理 。 


m= 一 0.1 == ¢=0.15 =m C=0.2 


(b) 阻尼 比 对 稳定 性 的 影响 


m= [og—(0.4 =m Ls=(0.5 mm Lo=0.67 


(d) 电 网 强度 对 稳定 性 的 影响 
6 系统 参数 对 稳定 性 的 影响 


Fig.6 Influence of different parameters on stability 
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mm fp =15Hz mm fp =25H7 = ee (一 0.2 ma (一 0.1 mm cosp=] mm cosp=0.6 mom /,=0.4 mmm L,=0.67 


1.005 . 1.002 1.002 1.2 
= = = = 
& 1 | = 1 &1 
x x x x 
0.995 0.998 0.998 . 0.8 
1 2 3 1 15 2 1 15 2 1 1.5 2 
时 间 /s 时 间 /s 时 间 /s 时 间 /s 


7 参数 变化 后 的 电压 响应 曲线 


Fig. 7 The voltage responses with various parameters 


42 ”和 鲁 棒 性 分 析 与 验证 

在 Matlab/Simulink 中 搭建 图 1 所 示 的 VSC 并 
网 系统 ， 参 数 见 附录 B 案例 1。 首 先 利 用 相同 摄 动 
下 开 环 传递 函数 奈 奎 斯 特 曲线 的 变化 验证 3.2、3.3 
节 的 结论 。 假 设 线路 中 的 电感 具有 乘 性 不 确定 性 ， 
摄 动 相 对 大 小 为 5%， 则 摄 动 后 的 电感 导 纳 大 小 为 
L=L,(l+6e), 9¢€|0,360 | (35) 
图 8 为 线路 电感 Ls 遭受 相同 摄 动 下 三 种 开 环 
传递 函数 La Lp 和 Lo 的 奈 硅 斯 特 曲线 带 ， 当 VSC 
运行 功率 因数 为 0.9 (140=0.90, Lo= 一 0.44), 其 裕 度 
变化 量 如 图 8 (a) 所 示 。 进一步 , 改变 外 环 控制 模式 
为 定 直流 电压 控制 ，VSC 运行 功率 因数 为 08 
Cao=0.80, I0=— 0.60), 其 裕 度 变化 量 如 图 8 (b) 所 
示 。 上 述 结果 表明 ， 在 受到 不 确定 摄 动 时 ， 广 义 阻 
抗 判 据 所 分 析 的 开 环 传递 函数 Lo 在 处 理 PLL 参与 

度 高 的 振荡 问题 时 鲁 棒 性 更 好 。 


mw ly eee Lp oe LO 


5 


© 


虚 
轴 


10 -12 -1 -08 -0.6 -0.4 
实 轴 


b) 直流 电压 外 环 控制 ， 功 率 因数 为 0.8 
8 考虑 摄 动 后 的 奈 奎 斯 特 曲 线 带 


Fig.8 Nyquist curve band under perturbation 
接 下 来 , 基于 三 种 开 环 传递 函数 整定 PLL 参数 ， 
说 明 在 相同 矢量 裕 度 下 VSC 并 网 系统 的 鲁 棒 性 存 
在 差异 。 设 置 VSC 运行 功率 因数 为 0.8 (1ao=0.80， 
17ao= 一 0.60)， 其 余 参 数 见 附录 B 的 案例 2， 系 统 的 
主导 振荡 模式 为 si= -1.24+181.7, 其 阻尼 比 为 0.0152, 


利用 La Lp 和 Lo 得 到 的 矢量 裕 度 分 别 为 0.389, 

0.166 和 0.145。 根 据 控制 理论 的 要 求 ， 矢 量 裕 度 应 
该 不 小 于 0.5P9]， 可 以 通过 整定 PLL 的 参数 来 提升 
系统 的 矢量 裕 度 ， 如 表 4 所 示 。 参 数 整定 后 ， 三 种 
判 据 对 应 开 环 传递 函数 的 矢量 裕 度 均 为 0.5, 在 tls 
F, VSC 端 电压 受到 0.01pu 扰动 ，VSC 输出 的 d 
由 和 9 轴 电 流 时 域 波形 如 图 9 (a)、(b) 所 示 ， 可 以 
发 现 均 能 够 快速 收敛 , 说 明 三 种 PLL 参数 下 系统 均 
有 一 定 的 稳定 裕 度 。 在 拓 3s 时 ,线路 电感 变 为 0.5Spu， 
通过 对 比 图 9 中 VSC 输出 的 电压 和 电流 波形 ， 发 
现 基 于 广义 阻抗 判 据 整 定 的 VSC 系统 鲁 棒 性 更 好 。 

表 4 基于 不 同 阻抗 判 据 的 PLL 参数 整定 


Tab.4 PLL parameters setting based on different criteria 


TH 


Ae} 


广义 阻抗 判 据 | dq 坐标 阻抗 判 据 序 阻抗 判 据 
0.5 
-6.81+j83.0 
mm lo 
-0.3 
-0.6 
~ 
-0.9 
0.6 1.2 
0 


2 3 4 0 1 2 3 
时 间 /s 时 间 /s 
(a) 输出 4 轴 电 流 (b) 输出 9 轴 电 流 
9 电磁 仿真 中 系统 振荡 波形 
Fig.9 Waveform in EMT simulation 


5 结论 


不 同 阻 抗 法 在 判断 VSC 并 网 系统 的 稳定 性 方 
面 等 价 ， 但 它们 的 标 称 性 和 和 鲁 棒 性 不 同 ， 这 导致 不 
同方 法 在 解决 相同 问题 时 的 效果 存在 差异 。 阻 抗 判 
据 的 高 鲁 棒 性 有 利于 VSC 并 网 系统 的 控制 环节 设 
计 和 参数 整定 , 广义 阻抗 判 据 在 PLL 参与 度 高 的 次 
/ 超 同步 振荡 问题 中 鲁 棒 性 较 好 , 序 阻抗 判 据 在 内 环 
和 线路 LC 参与 度 高 的 中 高 频 振荡 问题 中 鲁 棒 性 较 
好 。 针 对 PLL 主导 的 次 / 超 同步 振荡 问题 ， 当 功率 
因数 为 1， 广义 阻抗 和 dq 阻抗 判 据 的 标 称 性 等 价 ， 
然而 当 考 虑 VSC 多 象限 运行 场景 ， 不 同 阻抗 判 据 
的 标 称 性 存在 差异 ， 广 义 阻抗 判 据 所 分 析 的 开 环 传 
递 函 数 不 含 RHP 极点 ，dg 阻抗 判 据 和 序 阻抗 判 据 
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则 难以 避免 此 问题 

阻抗 法 本 质 是 一 种 特殊 的 频 域 法 ， 其 适应 性 和 
稳定 机 理 紧密 相关 ， 甚 至 很 多 场合 下 不 是 最 佳 的 频 
域 方 法 ， 故 寻找 更 适合 的 频 域 方法 以 及 这 些 方法 与 
问题 的 匹配 规律 是 未 来 需要 回答 的 问题 。 类 似 地 ， 
上 述 匹 配 性 的 问题 在 多 机 系统 中 也 存在 , 如 众多 短 
路 比 的 适应 性 问题 ， 这 也 是 未 来 需要 开展 的 工作 。 
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附录 A 
线性 化 文献 [25] 中 变 流 器 的 动态 方程 式 (1)-(9) 并 转换 到 以 同步 转速 旋转 的 电网 坐标 系 可 以 得 到 全 局 dg 


坐标 系 下 的 变 流 器 阻抗 模型 (为 表示 方便 ， 下 文 用 dg 表示 全 局 坐标 系 ， 电 流 正方 向 为 流出 变 流 器 ) 
Gyr lS) + Gi(S)G xo (San (Gou (9) + Gae (8)) GLa 
M __| GOGN, 1+ G,(8)Ga (SW, a (Al) 
Al, GCC)Cu(G) (I= Vo Gus. (8)) Gv) = (Gor. (9) + Gre (9))G (Ma [LAV 
1+G,(s)Gap(s WV, 1+G,(s)Gap(s WV, 


式 中 各 元 素 表 达 式 及 交流 电网 阻抗 模型 推导 详 见 文献 [25]。 


由 于 在 本 文 所 研究 的 频段 内 ，VoGpri(s)Gvr(s) 的 幅 值 很 小 ， 如 图 Al(a) 所 示 ,， 且 忽略 该 元 素 前 后 VSC 阻 
Hu REA Bode 图 基本 重合 ， 如 图 A1(b) 所 示 。 因 此 ，W0GpLL(s)Gvr(s) 可 近似 忽略 ， 上 式 可 以 近似 为 
| Gul) +G Gao (Cu +CD) Go | 
Als] | GOC HAO ar 
网 | CC (sMo Gy; (s) - (Gou (s) + Grr (s))G, (sao K a 
1+G,(s)Gap(sM% 1+ G,(s)Gap(s)V% 


式 中 : o, 1q0 分 别 对 应 d HH, q 轴 电 流量 的 稳 态 值 ， 在 单位 功率 因数 下 Tao=1.0 pu, 450 pus 
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Fig. A1 Bode diagram of the VSC q-axis transfer function 


附录 B 
RBI 变 流 器 并 网 系统 参数 收 稿 日 期 
Tab.B1 Parameters of the grid-connected converter 作者 简介 : 


参数 数 信 
E 
系统 容量 基 值 Sw/ KVA 1500 
直流 电压 基 值 Ubao V 1100 h 
直流 电容 Ca pu 0.038 N 


艾 力 西 尔 : 亚 尔 买 买 提 (1996), 男 , 硕 
研究 生 , 研究 方向 为 新 能 源 电力 系统 稳 
` 定性 分 析 和 控制 ，alxr@zju.edu.cn; 


* 通 信 作 者 : 辛 焕 海 (1981)， 男 ， 教 


滤波 电感 电容 Ze Cy pu 0.05,0.05 0.05,0.05 HDR WR 授 ， 博 士 生 导 师 ， 主 要 研究 方向 为 新 能 源 
aap 定性 Es 
电流 环 PI 0.8,18 0.6,8 买 买 提 电力 系统 稳定 性 分 析 与 控制 ， 
电压 前 馈 Te 0.001 0.002 a 
xinhh@zju.edu.cn. 


DEAT, 230 
锁 相 环 PI 18,8100 12,7200 


TE: 本 文 参考 文献 [27] 的 案例 参数 
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There are many impedance-based methods to an- 
alyze the oscillation of grid-connected converter sys- 
tem (GCS) via dq frame, sequence frame and polar 
frame, which can be transformed into each other and 
are equivalent in judging the stability of the system. 
However, the discrepancy and adaptation of these 
methods has not been discussed theoretically. In this 
paper, several impedance criteria of GCS were com- 
pared in terms of adaptation for analyzing the sub-syn- 
chronous oscillation(SSO) caused by the phase-locked 
loop (PLL), and the existence of open-loop transfer 
function(OLTF) right half plane(RHP) poles and its 
sensitivity to uncertain parameters are analyzed from 
the nominal performance and robustness aspects of the 
impedance criteria. 

A common method is to eliminate some variables 
via mathematical transformation. Then the complex 
multi-input multi-output(MIMO) system is reduced to 
an equivalent single-input single-output(SISO) system. 
Based on this, the stability analysis and control design 
could be carried out using methods such as Nyquist sta- 
bility criterion. To effectively guide the control design, 
the OLTF should satisfy the nominal performance and 
(1) without RHP poles. (2) the stability 
margin won’t have completely different result with un- 
certainty perturbations. 

The signal flow diagram of the equivalent SISO sys- 
tem based on impedance model in dg frame, sequence 
frame and polar frame when the outer loop and the volt- 
age feedforward filter dynamics were ignored is shown 
in Fig.1, which shows that the part dynamics of the PLL 
are difficult to separate from the network dynamics or 
power factor in L,(dq frame) and Lp(sequence frame), 
as depicted in Fig.1(a) and (b). This structural feature 
may lead to the existence of RHP poles and sensitivity 
to uncertain parameters in OLTF. In contrast, It can be 
seen that the PLL closed-form expression are cascaded 
with the the rest of the OLTF in Lo(polar frame), as 
shown in Fig.1(c). 


robustness: 


(a) Signal flow diagram of the La 


让 av, 


(c) Signal flow diagram of the Lo 

Fig.1 Diagram of the equivalent SISO system 
Further, the nominal performance and robustness 
of different impedance criteria is analyzed in the case 
of the SSO caused by the PLL considering converter 
four quadrant operation scenario. In nominal perfor- 
mance, it is found that the stability of the system can be 
accurately judged based on the information of RHP 
poles and Nyquist curves corresponding to the three 
impedance criteria, as shown in Fig.2(a). However, Za 
and Lp may have RHP poles in terms of this oscillation 
issue, which makes it difficult to quantify the system 
stability margin. In robustness, it is found that the sen- 
sitivity of the Leis smaller than La and Lp by comparing 
Nyquist curve band of open loop transfer function un- 

der the same perturbation, as shown in Fig.2(b). 
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(b) Nyquist curve band under perturbation 


Fig.2 the nominal performance and robustness analysis of 
the impedance criteria. 


